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All physical implementations of quantum bits 
(qubits), carrying the information and compu- 
tation in a putative quantum computer, have 
to meet the conflicting requirements of environ- 
mental decoupling while remaining manipulable 
through designed external signals. Proposals 
based on quantum optics naturally emphasize the 
aspect of optimal isolation while those fol- 

lowing the solid state route exploit the variabil- 
ity and scalability of modern nanoscale fabrica- 
tion techniques ||^^|. Recently, various designs 
using superconducting structures have been suc- 
cessfully tested for quantum coherent operation 
l^^lj], however, the ultimate goal of reaching co- 
herent evolution over thousands of elementary 
operations remains a formidable task. Protect- 
ing qubits from decoherence by exploiting topo- 
logical stability, a qualitatively new proposal due 



to Kitaev |12|, holds the promise for long deco- 
herence times, but its practical physical imple- 
mentation has remained unclear so far. Here, we 
show how strongly correlated systems developing 
an isolated two-fold degenerate quantum dimer 
liquid groundstate can be used in the construction 
of topologically stable qubits and discuss their im- 
plementation using Josephson junction arrays. 

Any quantum computer has to incorporate some fault 
tolerance as we cannot hope to eliminate all the various 
sources of decoherence. Amazing progress has been made 
in the development of quantum error correction schemes 
|p3| which are based on redundant multi-qubit encoding 
of the quantum data combined with error detection- and 
recovery steps through appropriate manipulation of the 
data. Error correction schemes are generic (and hence are 
applicable to any hardware implementation) , but require 
repeated active interference with the computer during 
run-time; the delocalization of the data, often in a hi- 
erarchical structure, boosts the system size by a factor 
10^ to 10"^. Delocalization of the quantum information 
is also at the heart of topological quantum computing 
| p2t , however, the stabilization against decoherence is en- 
tirely deferred to the hardware level (hence it is tied to 
the specific implementation) and is achieved passively. 
In searching for a physical implementation of topological 
qubits one strives for an extended (many body) quan- 
tum system where the Hilbert space of quantum states 
decomposes into mutually orthogonal sectors, each sec- 
tor remaining isolated under the action of local pertur- 



bations. Choosing the two qubit states from ground- 
states in different sectors protects these states from un- 
wanted mixing through noise; protection from leakage 
within the sector has to be secured through a gapped 
excitation spectrum. As no local operator can interfere 
with these states, global operators must be found (and 
implemented) allowing for the manipulation of the qubit 
state. 

A promising candidate fulfilling the above require- 
ments is the quantum dimer system fl^-jl^ : recent quan- 
tum Monte Carlo simulations of the dimer model on a 
triangular lattice provide evidence for a gapped liquid 
groundstate (see for a discussion of similar ex- 

otic groundstates in spin models) and we will discuss its 
topological robustness below. The physical implementa- 
tion of such a dimer system can be realized with the help 
of quantum Josephson junction arrays. Such arrays allow 
for a proper 'material design' as the energy scales govern- 
ing the charge and phase degrees of freedom can be tuned 
through junctions with appropriate tunnel barriers and 
capacitances. Shaped into the form of a ring, the Joseph- 
son junction array emulating a dimer system allows for 
the construction of topologically protected qubits. 

The simplest dimer model is defined on a square lattice; 
allowed configurations are coverings with dimers connect- 
ing neighboring vertices satisfying the constraint that ev- 
ery vertex belongs to only one dimer, see Fig. |^. Attribut- 
ing equal energies to these states, the classical model is 
characterized by an extensive entropy. The quantum dy- 
namics is introduced by the hopping Hamiltonian 
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rotating parallel dimers on appropriate plaquettes □ . As 
a result, the degeneracy of the classical model is lifted, 
with 'columnar' and 'staggered' (e.g., brickwall) phases 
(see Fig. ^ assuming minimal and maximal energies, re- 
spectively. In order to avoid such crystalline ordering we 
frustrate parallel dimers with the interaction term 



i/„ = «^(|ii)(ii| + |z)(z|). 
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At t ~ V (the Rokhsar-Kivelson point) a zero energy 
dimer liquid groundstate is formed involving all allowed 
configurations with equal amplitudes ]l5[ | . For the square 
lattice this liquid phase is restricted to the point t = v 
and tuning v away from t the liquid loses the competi- 
tion with various crystalline phases pS]. Here, we are 
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interested in systems developing a stable liquid over an 
extended parameter range. Such a stabilization may be 
achieved by the term 

H, = J2[-t{\_-) (//| + I, i> (^1 + h.c.) + Mk) (/| ] (3) 

□ 

transforming horizontally d-")) and vertically (|| ')) 
shifted dimers into diagonal ones |^) at the energy cost 
2fi (only one family of diagonals is chosen). Tuning the 
chemical potential fi from cx) to 0, diagonal dimer flips 
proliferate and we can study the evolution of the dimer 
model as we go from a square to a triangular lattice [El| . 




FIG. 1. Dimers on a square lattice forming liquid (main 
panel), brickwall, and columnar configurations. Periodic 
boundary conditions are applied along the 2;-direction. The 
sublattice A is marked by solid bullets; dimers are attributed 
a direction pointing away from bullets. The topological in- 
variant counts the number of dimers crossing the reference 
line 7, with dimers pointing to the right (left) counting with 
positive (negative) signs (here, </!> = 1). Right: staggered and 
columnar configurations in the triangular lattice. 

The above type of quantum dimer models has been ac- 
tively discussed in the context of short-range resonating- 
valence-bond (RVB) models for high temperature 
superconductivity In order to gain more insight 

into its physical properties we exploit a mapping to an 
'electrodynamic' analogue po| , ^ , p3| : We decorate the bi- 
partite square lattice (sublattices A and B) with static 
charges pn = ±1, placing all positive charges on the 
sublattice A, see Fig. |l|. A dimer covering {d^j j^} of 
bonds of this 'ionic crystal' is mapped to a dis- 

tribution of electric dipoles {— p^j^-^}; the electric field 
e = —{p{ij))sA resulting from averaging over small ar- 
eas SA then satisfies Gauss' law Ve = {p)5A = 0. The 
simplest dynamics preserving this constraint is given by 
Faradays law or, more precisely, (compact) quantum elec- 
trodynamics. Here, we are interested in the large-scale 
(global) aspects of the model and hence we will keep 
the discussion on a phenomenological level; the physical 
properties of the dimer liquid phase will be described in 
terms of the electric field e and the dielectric constant e. 

While the order characterizing the crystal phases is ob- 
vious, the topological order present in the dimer liquid 
phase is more subtle [p"4|-p^: we consider a cylindrical 
geometry with periodic boundary conditions along the 



X-axis, leaving the boundaries along the y-direction free 
(an additional periodic boundary condition along the y- 
direction defines a torus, see later). In the absence of 
diagonal dimers (the limit — s- oo in (^) the sum (j) 
counting directed dimers cut by the reference line 7 par- 
allel to the y-axis defines a topological order parame- 
ter (dimers pointing into opposite directions are counted 
with opposite sign) : the action of the Hamiltonian H = 
Ht + Hy leaves (j) invariant and the Hilbert space splits 
into topological sectors Ti^ characterized by the integers 
(j). Within our electrodynamic analogue this topological 
invariant is given by the flux cf) = J dy ee^ of the elec- 
tric induction across 7 and each sector Ti.^ is character- 
ized by an average polarization Cx = (t>/sLy producing a 
groundstate energy E = ee'^LxLy/2 — (j)^Lx/2eLy = E^. 
At finite values of the chemical potential p the Hamilto- 
nian (^) introduces mobile defects with charge ±2 |^ , p4| . 
These charges mix states with electric fiux </) ± 2: creat- 
ing two defects of charge ±2, carrying one defect around 
the cylinder, and subsequent annihilation with its part- 
ner changes the flux cj) by 2; as a result, the infinite set 
of sectors TY^ collapses to the even and odd sectors Tie,o- 

The Hamiltonian projected onto the space of ground 
states within topological sectors Ti^ takes the form 
^-c = E0(^0l<^>(0l + J{\<t>){^ + 2| + h.c.)). For large 
values of the chemical potential fj, ^ t the defect pairs 
are only virtual and the amplitude J mixing states with 
different (j) is exponentially small in the circumference L^, 
J (X {t/p)^". 

As ^ drops below an energy of order i, defects pro- 
liferate and undergo Bose condensation; the concomi- 
tant long range order in the condensate phase Lp im- 
plies large charge fluctuations and thus the electric fleld 
flux across 7 fluctuates strongly. In this phase the mix- 
ing between sectors is large with (t^^) cx Ly: assum- 
ing a superfluid density ps^ the hydrodynamic descrip- 
tion of the charged superfluid |^3| is given by the Eu- 
clidean action (in Fourier space; p denotes the conju- 
gate to ip) 5b = X]i^,K['^*°-'-',-K</'i^,K + Ps^^|'/'t^,Kp/2-|- 

|Pi^,Kp/2ei^^] , from which we derive a density correlator 
{pp) — PsK"^ / [Ps/^ + w^). The relation eV • e = p allows 
us to find the equal time correlator of the electric induc- 
tion (£e^(R,t)ee^(R',i)) = (5^j.(52(R-R')^p;?/2, result- 
ing in the mean electric flux (0^) = pssLyjl crossing 
the line 7; the local correlations in the electric field are 
a consequence of the dispersion-free plasmon spectrum. 

The new groundstates |e, o) oc ^ a^^ ol'/'e,o) involve 
contributions from all sectors 7ie,o, with 0e — 2k and 
(j)o = 2k 1 running over even and odd integers; their 
physical properties are conveniently characterized by the 
distribution function P{<j)) — ja^P oc exp(— 0^/2cr), with 
the Gaussian form appropriate for a process involving a 
sum of statistically independent elements. The width a 
is obtained through a comparison with the second mo- 
ment ((/)^), a — ^ps£Ly/2. Expectation values in even 
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and odd sectors 7ie,o are given by the sums over even and 
odd (j> with weight P{4>)- As a result, differences between 
expectation values in 7ie,o are exponentially small in Ly, 
e.g., the average squares of electric fluxes differ by 

(02),_(c/,2)„^exp(-7r2^L^/4). (4) 

Hence, the presence of diagonal dimers provides us with 
only two topological sectors for a system defined on a 
cylinder (and 4 sectors on a torus). The correspond- 
ing groundstates are degenerate in the limit L^^y oo. 
Finite-size effects as well as local disorder lead to an ex- 
ponentially weak energy splitting (cf. (^). Furthermore, 
the topological property of the system inhibits any mix- 
ing of the two groundstates. In the thermodynamic limit 
the collapse of infinitely many to two topological sectors 
involves a quantum phase transition in the parameter /i. 
For /i = the above model corresponds to a symmet- 
ric triangular lattice (upon including the v-term) and we 
expect favorable conditions for the construction of pro- 
tected qubits. 

Indeed, this expectation is supported by results ob- 
tained from numerical studies: recent Monte Carlo sim- 
ulations by Moessner and Sondhi [jl^ on large systems 
with Lx^y = 36 exhibit short-range dimer-dimer correla- 
tions within a parameter region 2/3 < v/t < 1, indicative 
of a liquid state. Furthermore, the weak temperature de- 
pendence of this data within the temperature interval 
0.25i > T > O.QSt provides evidence for a gapped spec- 
trum (increasing v beyond t, a first-order transition takes 
the liquid into the staggered phase, while decreasing v 
below 2t/3 establishes some crystalline order). 

We have numerically diagonalized the t-v dimer Hamil- 
tonian H = Ht + Hy on a triangular lattice in order to 
estimate the gap protecting the liquid and verify the ab- 
sence of low lying edge states; in addition we have in- 
vestigated the robustness of the degenerate groundstates 
under local perturbations (requiring exact diagonaliza- 
tion limited to smaller systems). We have chosen both 
cylindrical and toroidal geometries, with system sizes go- 
ing up to Lx^y = 6; although subject to finite size effects 
we expect that our approach describes well the liquid 
phase where the short-ranged dimer correlations (of or- 
der of one lattice constant ||l^) reduce the influence of 
the boundaries. Our exact diagonalization study con- 
firms the presence of a gap of order A « O.lt, taking up 
large values on approaching the Rokhsar-Kivelson point 
t — V. Comparing the spectra for tori (no open bound- 
aries) with those of cylinders (where the two open bound- 
aries could accommodate edge states, would they exist) 
we note that the gap persists, from which we infer the 
absence of low lying edge-states in the dimer liquid (this 
should be contrasted with the situation in the quantum 
Hall system). 

In order to test the susceptibility of the degener- 
ate groundstates to local perturbations we assign all 



links random chemical potentials Hd homogeneously dis- 
tributed over the interval [—d/2,d/2], d < t. In finite 
cylinders and tori, particular symmetries specific to some 
geometries L^.y produce an exact degeneracy of eigen- 
states from different topological classes; for the 6*5 cylin- 
der and torus used in our study the symmetry analysis 
predicts double-degeneracy. As shown in Fig. ^the disor- 
der induced groundstate splitting collapses dramati- 
cally as we go from the crystal {v > t) to the liquid phase 
{v < t) and slowly recovers as v is decreased further. We 
attribute the sharp drop at v » t to the first-order solid- 
liquid transition; the small disorder splitting a.t v < t 
testifies for the efficient protection of the dimer liquid 
groundstates from local perturbations. Decreasing v fur- 
ther below t we observe level crossings which we attribute 
to the appearance of intermediate crystalline phases |l7| ] 
before reaching the fluctuating columnar phase at large 
negative values of v. 
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FIG. 2. Splitting of the groundstate energies under the 
action of a disorder potential of strength d for a 6 * 5 torus 
and a 6*5 cylinder (inset; periodicity is along x with Lx = 6)- 
The disorder only weakly perturbs the dimer liquid ground- 
states at 11 < t (note the expanded view near v ^ t for the 
6*5 torus). The crossovers between various phases become 
clearly visible in the susceptibility to a disorder potential. 

Summarizing, the combination of Monte Carlo sim- 
ulations fr^ and diagonalization indicate that on ap- 
proaching the point v = t from below, the triangular 
dimer model develops an isolated degenerate dimer liq- 
uid groundstate free of edge states. Within the param- 
eter region ~ 0.8 < v/t < 1 we find a gap of order 
O.lt and groundstates which are robust under disorder 
to within 1% of the disorder potential; we expect a fur- 
ther increase of this robustness with system size, cf. (Q). 
Translating these findings to potential implementations 
of this two-state system, we find the latter to be quite 
tolerant with respect to (static) local variations in pa- 
rameters appearing in the fabrication process. Regarding 
decoherence, the amplitude mixing of the two ground- 
states is strongly suppressed, since it involves creation of 
topological defects (real or virtual) violating the dimer 
constraint. The decoherence in the relative phase of the 
two groundstates originates either from their adiabatic 
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splitting by an external low-frequency noise, or the cre- 
ation of non-topological excitations within each sector. 
The former is suppressed for the same reason that leads 
to the robustness with respect to local disorder, while the 
latter is inhibited by the gap A, which sets the require- 
ment T ^ A on the operation temperature. 

In the following we propose two types of Josephson 
junction arrays for the implementation of the above topo- 
logically protected quantum dimer liquid groundstates 
and their use for quantum computing. In the first ar- 
ray ( J JT) , the vertices of the triangular lattice are struc- 
tured into six Y-shaped superconducting islands with two 
ends joining into the hexagonal vertex and the third end 
linking to the neighboring hexagon, see Fig. |^. The 5 pa- 
rameters, the capacitance Cy of the island to the ground- 
plate, the capacitance Ch and Josephson current Ih of the 
hexagonal vertex junction, and the capacitance C; and 
Josephson current /; of the link junction, are chosen to 
emulate the triangular quantum dimer model. All capac- 
itances and currents define corresponding energies E'^ — 
e^/2C and i?' = ^oI/2ttc, where $o = hc/2e denotes the 
flux quantum (e.g., = {2ef/2Ch, E^ = <Poh/2nc). 
We first find Cy, Ch, and /; defining the classical dimer 
states: We choose a large capacitance Ch in order to 
join the islands electrically into one hexagonal vertex. A 
small capacitance Cy defines the large charging energy 
-E'hcx ~ -E-y/S of the hexagonal vertex, the basic energy 
scale of the array. We introduce 'charge frustration' of 
the vertices by biasing the array with a global electric 
gate such as to equalize the energies (to an accuracy bet- 
ter than £"/) of two states differing by one Cooper pair 
(this defines a 'half- filled' array); the large charging en- 
ergy ii^hox lifts other charge states to high energies. The 
Cooper pairs lower their energy via tunneling (involving 
the coupling //) through the link junction joining two ver- 
tices — the 'bonding' state of such a Cooper pair defines 
the dimer state. With only half a Cooper-pair available 
per hexagonal vertex, each vertex is involved in the for- 
mation of one and only one of these 'valence bonds'. This 
defines the 'classical' configurations of the dimer model; 
the corresponding Hilbert space of dimer states is pro- 
tected by the energy scale Ef . 

The dimer dynamics involves the vertex junction and 
proceeds via localization of one dimer to the vertex (with 
an energy cost E^) and subsequent hops (of two parallel 
dimers/Cooper pairs over junctions with energy see 
Fig. |) to the new vertex islands. The hopping amplitude 
is of order t ^ E'^^ /E^; explicit calculation gives the re- 
sult t — {%/\&)E-^ /El. The electrostatic interaction be- 
tween dimers depends on the choice of the capacitance 
matrix and its dependence on the dimer configuration is 
non-trivial. First, we compare the energies of the stag- 
gered and columnar configurations and find the interac- 
tion energy v between parallel dimers; in the limit C; ^ 
C^^Ch, V = Elid/Chm + CYlCh){l + 3Cy/C,,)]-2. 
Second, we have checked that the electrostatic energies 



of the liquid phase dimer configurations indeed scale (to 
within ~ 5% accuracy) with the number of parallel dimer 
pairs. The condition C; ^ Cy guarantees a short range 
interaction between dimers. 
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FIG. 3. Josephson junctions arrays forming triangular 
(JJT) and Kagome (JJK) lattices and emulating dimer sys- 
tems. JJT: The large capacitance Ch joins hexagonal is- 
lands into electric units which are frustrated by a global 
gate to accept 1/2 Cooper pair. Cooper pairs resonat- 
ing between hexagonal islands form dimers (shaded ellipses); 
dimer flips due to Ht axe indicated by arrows. JJK: Cooper 
pairs on X-shaped islands are capacitively coupled through 
star-shaped islands. The gate is tuned to allow for 1/2 Cooper 
pair per hexagon. Dimers are pairs residing on X-shaped is- 
lands and polarizing adjacent star-shaped islands. 

The backbone of the second array (JJK) is made from 
X-shaped islands arranged into a Kagome lattice with 
capacitive and Josephson couplings Ef^ and E'l^, see Fig. 
^. A second triangular lattice of star-shaped islands in- 
troduced into the hexagons of the Kagome lattice is only 
capacitively coupled (with energy E'f) to the islands of 
the Kagome lattice; Cx and C* denote the capacitances 
of the islands to the ground. The Kagome sites are bi- 
ased with a global external gate to accept one additional 
Cooper pair per hexagon at equal energy, resulting in a 
filling with 1/2 Cooper pair per hexagon. The similar- 
ity between the Kagome lattice and the triangular array 
becomes obvious when joining pairs of Y-shaped islands 
(corresponding to the limit C/ ^ oo and // oo) and 
contracting the link to obtain the X-shaped island — the 
previous vertex junctions now play the role of the junc- 
tions on the Kagome lattice, hence the same index '/i' 
has been chosen for the junctions on the hexagons; note 
that the present array involves only one type of Josephson 
junctions. By analogy, dimers now correspond to Cooper 
pairs localized onto the X-islands of the Kagome lattice. 
As dimers should not touch one another, no two bosons 
are allowed on the same hexagon: We choose small ca- 
pacitances Ch in order to isolate hexagons from one an- 
other (note the difference to the JJT array) and large 
capacitances C^ to join the six X-shaped islands forming 
the hexagon into one electrostatic unit via their strong 
coupling to the central star-shaped island. The capaci- 
tances Cx, C*, and Ci then define the charging energy 
required to put two Cooper pairs on the same hexagon, 
Ehcx ~ {Ci/Cx)'^E^, which is our basic energy scale (we 
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assume Ci < C* , Cx in order to guarantee good screening 
on large distances). 

The motion of the Cooper pairs (the dimer dynam- 
ics) involves the hopping amplitude i?^ and the charging 
energy E'hex of the virtual state with two Cooper pairs 
on the same hexagon, t ~ E'l^ / Ehox', note that we re- 
quire Ef^ <C i?hex not to perturb the proper frustration of 
the array. The interaction energy v again derives from a 
comparison between the staggered and columnar states 
and has to meet the conflicting requirements of increas- 
ing the energy of parallel dimers (next-nearest neighbor 
interaction) while leaving non-parallel dimer pairs (next- 
next-nearest neighbor interaction) approximately unper- 
turbed. In order to produce a positive v we need a finite 
coupling Cfi, which in turn reduces the basic energy E-^cx- 
We find that the set 5Ci = = Cx = 20Ck produces 
a (renormalized) protective energy E^^^ « O.OOSE'J ~ 
0.2i?hex and a repulsive energy v « 0.1i?JJg^. Choosing 
Ef^ « O.SS^g^ we can properly satisfy the condition t fa v. 

We construct a topologically protected qubit from the 
two-level system defined by the two ground states of a 
quantum dimer system with cylindrical boundary con- 
ditions. We make use of the above Josephson junction 
arrays and a ring geometry; the two qubit states |e) and 
|o) then are distinct through the parity of the dimer count 
along the line 7 joining the inner and outer boundaries 
of the array, see Fig. ^ In order to manipulate the qubit 
we have to implement the qubit Hamiltonian 

^^qubit = hxCTx + hz(Tz, (5) 

where Ux and are Pauli matrizes and h^, hz are 
the (manipulable) parameters producing the amplitude 
(a) and phase (x) mixing in the qubit state \a,x) = 
[ |e) -|- a exp(zx)|o) ]/VT + c?. The implementation of hx 
requires a controlled mixing of the protected dimer states, 
implying a reduction of the ground state's topological 
protection. In the JJT array this is achieved by break- 
ing one dimer bond and creation of a virtual particle- 
hole excitation where one Cooper pair retreats to one 
hexagon (particle), leaving the partner hexagon empty 
(hole). This virtual excitation costs the energy E^ of 
that particular bond and takes the system (virtually) out 
of the protected dimer space. While the particle remains 
pinned to the weak junction, the hole is taken around the 
inner boundary through appropriate dimer flips and is 
subsequently recombined with the particle. This process 
results in a mixing amplitude hx ^ El^{El^/ E^Y'^ , where 
M denotes the number of links on the inner boundary 
(this estimate applies to the 'optimal' dimer configura- 
tion shown in Fig. Hence introducing one switchable 
link junction near the inner boundary allows us to tune 
the parameter E'l and change the mixing amplitude hx 
by many orders of magnitude. The variant for the JJK 
array involves a virtual excitation with two Cooper pairs 
on one hexagon. 




FIG. 4. a The topologically protected qubit involves the 
two degenerate groundstates of the quantum dimer liquid 
state realized within a triangular Josephson junction array. 
The qubit 's phase is manipulated through the gated super- 
conducting strip (GSS) attracting dimers along 7. The am- 
plitude shifter mixing the groundstates is realized with the 
help of a tunable Josephson junction placed at the inner qubit 
boundary. Dimer-flips, like a row of falling dominos, take the 
hole around the inner boundary, b Array of topologically pro- 
tected qubits embedded in a Josephson junction array emulat- 
ing a dimer system. The entanglement of two qubits proceeds 
through coupling over a gated superconducting strip. 

The 'phase shifter' can be implemented through a 
gated superconducting strip, capacitively coupled to the 
array and attracting dimers onto the reference line 7. 
Upon biasing the strip, the energies of the two ground- 
states are shifted with respect to one another and the 
qubit 's phase x is modified accordingly. The implemen- 
tation requires careful adiabatic switching, with the am- 
plitude u (the energy acting on one dimer) and the du- 
ration r of the manipulation limited by the constraints 
u ^ t El and t > h/A in order to avoid excitations 
within the dimer liquid. Furthermore, we have to break 
up the strip during idle time, as the fully connected strip 
represents a global operator (otherwise, electric fluctua- 
tions fed to the strip would decohere the system). Thus 
the strip has to be constructed from isolated supercon- 
ducting islands which are to be connected only during the 
time where the phase shifter becomes operational. Fast 
switches accomplishing such a connection can be built 
from superconducting Cooper pair transistors p5|] . 

The qubit 's manipulation involve strong modifications 
of the system: the amplitude shifter requires leaving the 
protected Hilbert space through virtual breaking of a 
dimer, while the phase shifter requires introduction of 
a global operator. In turn, these processes are strongly 
inhibited during idle time and give the qubit its robust- 
ness: amplitude mixing is exponentially small in Lx while 
phase drifting is exponentially small in Ly, cf. (^. 

Constructing a register of qubits is surprisingly simple: 
the same basic Josephson junction array can accommo- 
date K qubits in a geometry with K holes, see Fig. |^. 
Individual qubit operations are carried out as described 
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above. The implementation of a non-trivial two-qubit 
operation again involves a superconducting strip: bias- 
ing the strip connecting two qubits lifts the energy of the 
states |eo) and |oe) with respect to the states |ee) and 
\oo) and combining this two-qubit 'phase shifter' with 
a suitable set of single qubit operations allows for the 
construction of the controUed-NOT operation |^ . 

In a first step towards the implementation of topolog- 
ically protected qubits one may wish to test the proper 
functionality of the quantum Josephson junction array. 
Such a test is provided by a measurement of the mag- 
netic susceptibility x of the ring-shaped structure: ap- 
plying a magnetic field, the mixing energy hx picks up a 
factor cos(27r<l>/<i>o), where <I>o = hc/e is the normal state 
flux quantum. This factor originates from the Aharonov- 
Bohm phase picked up by the virtual charged excitations 
encircling the hole. In the absence of disorder the ring's 
energy changes with flux, i?i.ing($) = cos(27r$/$o)|, 
resulting in a suceptibility x (?|,£'iing with period 
$0/2 and a sharp feature at $ = <i>o/4. The split- 
ting Arf in the groundstate energies (due to disorder or 
finite-size effects) smears this sharp feature over a region 
(5<I)/$Q ^ Ad//ix; hence a measurement of x(^) allows 
to check on the mixing amplitude hx and on the level 
splitting. 

A particular challenge faced by qubit implementations 
based on Josephson junctions operating in the charge 
limit is the presence of stray charges [DJ^; the robust- 
ness in our device efficiently suppresses fluctuations in 
the gate potential once they drop below t. The problem 
of stray charges is avoided in designs using the dual phase 
variable and work in this direction is in progress. 
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